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Abstract 

We have calculated the chemical trend of magnetic exchange parameters (J^d, Na^ and Nf3) 
of Zn-based II- VI semiconductors ZnA (A=0, S, Se, and Te) doped with Co or Mn. We show 
that a proper treatment of electron correlations by the LSDA+fJ method leads to good agreement 
between experimental and theoretical values of the nearest-neighbor exchange coupling J^^ between 
localized 'id spins in contrast to the LSDA method. The exchange couplings between localized spins 
and doped electrons in the conduction band Na are in good agreement with experiment as well. 
But the values for Nf3 (coupling to doped holes in the valence band) indicate a cross-over from 
weak coupling (for A=Te and Se) to strong coupling (for A=0) and a localized hole state in 
ZnO:Mn. That hole localization explains the apparent discrepancy between photoemission and 
magneto-optical data for ZnO:Mn. 

PACS numbers: 75.50.Pp,71.23.An,71.55.Gs 
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I. INTRODUCTION 



After the seminal discovery of ferromagnetism in GaAs:Mn- with a critical temperature 
Tc as high as 110 K there is worldwide a renewed interest in diluted magnetic semiconductors 
(DMS). Recently, the Curie temperature in GaAs:Mn could be pushed to values of about 
180 K by a careful control of the annealing conditions during the growth process.™ There 
is a great search activity to look for alternative materials, especially in the class of II- VI 
semiconductors (SC). Ferromagnetism (FM) in diluted II- VI SC is known for a long time 
with up to now low values, however.- They also serve as model materials since they allow 
to control the magnetic ions and the doped charge carriers independently. In such a way 
it was possible to demonstrate the carrier-induced mechanism of the ferromagnetic state in 
Pb-doped SnTe:Mii^ or in p-doped ZnTeiMn.- 

The DMS combine ferromagnetism with the conductivity properties of semiconductors. 
Therefore, they are ideal materials for applications in spintronics where not only the electron 
charge but also the spin of the charge carrier is used for information processing. For instance, 
they allow to resolve the conductivity mismatch problem which hinders a high polarizability 
of injected electrons in a ferromagnetic metal/semiconductor junction.- 

The ferromagnetism in the traditionally known DMS arises due to Zener's p-d exchange 
mechanism.'' The 3d transition metal impurities lead to localized spins Sj. Hole doping into 
the valence band (either by the 3d transition metals itself or by other acceptor impurities) 
provides charge carries whose spins interact with the 3d spins. This local p-d exchange 
coupling = N(3 leads to a parallel arrangements of the magnetic moments since a fer- 
romagnetic state allows a higher mobility of the doped holes. For a high doping level the 
material becomes more metallic and the mechanism changes to a RKKY-like interaction. 

From this argumentation follows immediately that the crucial parameter to increase Tc 
is the Jp^ coupling. Indeed, a simple theory of Zener's p-d exchange mechanism^ gives 
Tc oc {Jp^)^Xh where Xh is the hole doping level. It can be expected from general grounds 
that a decreasing anion-cation distance leads to an increase of the p-d tight-binding hopping 
parameter tpd, and consequently to an increase of J^^. That reasoning lead Dietl et to 
the proposal of room temperature ferromagnetism in Mn-doped ZnO or GaN, respectively, 
which created a tremendous activity and numerous reports on room temperature FM in 
II- VI DMS or similar materials.-ii^iii 
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However, there are serious doubts whether the reported room temperature ferromag- 
netism belongs really to the same class of ferromagnetism as that one observed in GaAs:Mn 
or ZnTe:Mn which is based on Zener's p-d exchange mechanism. For instance, in ZnO:Co 
ferromagnetism was reported in samples produced by laser ablation,— or by the sol-gel 
method,^ whereas other samples fabricated by precursor deposition, or molecular beam 
epitaxy (MBE)^'^ showed no signs of ferromagnetism and antiferromagnetic couplings be- 
tween nearest neighbor 3d spins. It is highly probable that the observed ferromagnetic 
effects in ZnO:Co are due to uncompensated spins at the surface of Co- rich antiferromag- 
netic nanoclusters.— 

The proposal of Dietl et af was based on simple model calculations and qualitative 
arguments. There is a real need for a parameter free ab-initio study of the relevant exchange 
parameters in II- VI semiconductors to put the expected chemical tendency on a firm basis. 
Such a calculation of the nearest neighbor couplings of local spins Jdd and the p-d exchange 
couplings and J^^ with valence and conduction bands, respectively, is presented here. 
We considered the series of Co- and Mn-doped ZnA with the anions A=Te, Se, S and O. 

To achieve our goal we had to solve two theoretical problems. First of all, the local spin 
density approximation (LSDA) is not sufficient. It leads to wrong predictions of FM in 
ZnO:Co even without additional hole doping,— to too large values of \Jdd\ for ZnO:Mn, and 
to the wrong (FM) sign of one of the two nearest neighbor exchange couplings in wurtzite 
ZnO:Co.^'' It was shown that this deficiency of LSDA can be repaired by taking into account 
the strong Coulomb correlation in the 3d shell by the LSDA+f/ method. To choose the U 
values we have to take into account that they decrease in the series from O to Te due to an 
increase of screening effects. The values of Jdd are very well known experimentally in this 
series. Therefore they can be used to check the chosen U values. We will show below that 
for reasonable values of U we obtain Jdd in good agreement with experimental results and 
we may explain the chemical tendency. 

The second theoretical problem concerns the p-d exchange coupling between the localized 
spins and the holes in the valence band J^^. This coupling leads to the giant Zeemann 
effect— and it is seen in our calculations by a band-offset AE"" between spin up and spin 
down of the valence band. For small values of J^^ (which means also small values of tpd) 
both splittings, the experimental and the theoretical one, are proportional to the magnetic 
impurity concentration x. In that weak coupling regime the p-d coupling can be simply 
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calculated by using the proportionality between splitting and x. We will show, however, 
that there are more and more deviations from AE^ oc a; if we go from ZnTe to ZnO. The 
exchange values obtained in that manner seem to depend on the concentration of magnetic 
impurities. We solve that problem by a fit to the Wigner-Seitz approach of Benoit a la 
Guillaume et al}^ Our results prove that we reach the strong coupling limit for ZnO. As we 
will show below, in that case the impurity potential is so strong that it can bind a hole for 
ZnOiMn, whereas ZnO:Co is close to the localized limit. 

Our ab-initio results strengthen the recent model calculations in Ref. |20|. That work was 



aimed to explain the tremendous difference between the experimental J^^ values obtained 
from photoemission and magneto-optics, especially in ZnO and GaN.^^ It was argumented 
that this difference arises due to state localization which is confirmed by our ab-inito calcu- 



lations below. But we also will show that our results for \N(3\ are much smaller than t 



hose 
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evaluated earlier from photoemission for ZnO:Mn (-2.7 eV (Ref. |2z) or -3.0 eV (Ref. 
and ZnO:Co (-3.4 eV (Ref. [2^) and which were used as model input parameters in Ref. 

The organization of our paper is as follows. After presenting the super-cell method in Sec. 
11, we discuss the nearest neighbor exchange coupling in Sec. III. That fixes the U values 
unambiguously. In Sec. IV we present our results for J'^^ = N(3 and J^^ = Na. Finally, in 
Sec. V we discuss the arguments in favor of a localized state in ZnO:Mn. 



II. SUPER-CELL CALCULATIONS 



We used super-cell calculations to determine the exchange couplings J^d, Na, and N(3. 
Since we are mainly interested in the chemical tendency within the II- VI series we restrict our 
study to the zinc-blende structure. All compounds of the series exist in that modification, 
even ZnO as epitaxial layer. To calculate Jdd we used super-cells of the form T2Zn6A8 with 
the transition metals T=Co or Mn and with the anions A=0, S, Se, and Te. In those 
super-cells the magnetic ions build chains. The exchange constants are then determined by 
comparing the total energies of ferro- and antiferromagnetic arrangements. We have checked 
that the infiuence of finite size effects is negligible (not larger than 6 per cent for Jdd) by 
performing some calculations with T2Zni4Ai6 super-cells. 

For we used super-cells with three different concentrations of magnetic ions, x = 1/4, 
X = 1/8, and x = 1/32, i.e. TZn3A4, TZnyAg and TZn3iA32. As will be explained below. 
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these numerical results have to be fitted with the Wigner-Seitz approach to obtain J^^. In 
all calculations we used the experimental lattice constants a = 6.101 A, 5.668 A, and 5.410 
A for ZnTe, ZnSe, and ZnS, respectively.— For ZnO we used a lattice constant a = 4.557 A 
which gives the same unit cell volume as the experimental value. (Bulk ZnO has a = 3.2427 
A and c = 5.1948 A in the wurtzite structure.)^ 

The super-cell calculations were performed using the full-potential local-orbital (FPLO) 
band structure scheme.— In the FPLO method (version FPL05) a minimum basis approach 
with optimized local orbitals is employed, which allows for accurate and efficient total en- 
ergy calculations. For the present calculations we used the following basis set: Zn,Co,Mn: 
3s3p:4s4p3(i, O: 2s2p;3d, S: 3s3p3d, Se: 4s4p3(i, and Te: dsdpM. The site-centered poten- 
tials and densities were expanded in spherical harmonic contributions up to Imax = 12. 

The exchange and correlation potential was treated in two different ways. First, the local 
spin-density approximation (LSDA) was used in the parametrization of Perdew and Wang.— 
However, as will be shown below in more detail, this approximation has severe deficiencies 
in the present case. The energetical positions of the Co(Mn) 3d states with respect to 
the valence band are incorrectly given in the LSDA calculation. They are expected to be 
much lower in energy and this correlation effect was taken into account by using the FPLO 
implementation of the LSDA+[/ method in the atomic limit scheme.— The convergence 
of the total energies with respect to the k-space integrations were checked for each of the 
super-cells independently. The calculations for each cell were first performed within the 
LSDA approximation using basis optimization. The LSDA+[/ calculations were then made 
starting from the LSDA optimized basis but with no basis optimization in the self-consistency 
cycle in order to obtain convergence. The Slater parameters and for Mn and Co in the 
LSDA+f/ calculations were chosen close to atomic values, namely = 7.4 eV and F^ = 4.6 
eV (corresponding to the Hund exchange coupling Jh = 0.86 eV) for Mn, and F^ = 7.9 eV 
and = 5.0 eV {Jh = 0.92 eV) for Co. The Slater parameter F° = U, however, is much 
more screened and its infiuence has been investigated more in detail (see below). 

III. D-D EXCHANGE COUPLINGS 

In this Section we are going to determine the exchange couplings between two localized 
magnetic ions. We are considering two nearest neighbor impurities, each carrying a local 
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FIG. 1: (Color online) Calculated exchange couplings Jdd for ZnA:Mn (from above to below: A=0 
(blue), S (red), Se (green), and Te (yellow) as a function of the Coulomb correlation U in the 3d 
shell. 
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FIG. 2: (Color online) Coulomb correlation U in the 2>d shell of Mn-impurities in II-VI-SC as func- 
tion of charge transfer 5Q. The values for CdTe, CdS, and ZnO are taken from Ref. I32 (constrained 
DFT-calculations, blue diamonds) and the other values (red squares) by linear interpolation cor- 
responding to the calculated charge transfer. 

spin Sj. Then, the Heisenberg Hamiltonian for a locahzed pair of spins is given by 

H = —2JddSiSj . (1) 

The corresponding total energies per magnetic ion for ferromagnetic (FM) and antiferromag- 
netic (AFM) arrangements of the two spins, EpM and Eafm, lead to the energy difference 
between the FM and AFM states: 

AE = ^^^^-^^^^^ = -^SriSr + 1) , (2) 
6 



TABLE I: Comparison of calculated and experimental values for the nearest neighbor exchange 
couplings Jdd for ZnA:Mn. 





exp. 






T A 




J,7(meV) 




f/(eV) 


Jrfd(meV) 


ZnO:Mn 


-2.09'' 


-2.18 


6 


-6.02 


ZnS:Mn 


-1.41'' ; -1.39'^ 


-1.39 


6 


-5.41 


ZnSe:Mn 


-l.OS"^ ; -1.06'= 


-1.19 


4 


-3.25 


ZnTe:Mn 


-0.75^ ; -0.76^ 


-0.65 


4 


-2.29 



'^Magnetization step method, Ref. 
''Magnetization step method, Ref. 
'^Inelastic neutron scattering, Ref. 
'^Magnetization step method, Ref. 
'^Magnetization step method, Ref. 
•^Inelastic neutron scattering, Ref. 



33, first neighbor in the (a,b) plane of the wurtzite structure. 
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TABLE II: Comparison of calculated and experimental values for the nearest neighbor exchange 
couplings Jdd for ZnA:Co. 





exp. 


LSDA+C/ 




LSDA 




J,^7(meV) 


Jdd(meV) 


;7(eV) 


Jdd(meV) 


ZnO:Co 


-2.0" 


-1.73 


6 


-1.39 


ZnS:Co 


-4.09** 


-4.13 


4 


-7.26 


ZnSe:Co 


-4.26'= 


-3.36 


4 


-6.26 


ZnTe:Co 


-3.27'= 


-3.32 


4 


-6.94 



'^Inelastic neutron scattering, Ref. 3§, first neighbor in the (a,b) plane of the wurtzite structure. 
^'Inelastic neutron scattering, Ref. 
"^Inelastic neutron scattering, Ref. |41 



where St is the total spin of two parallel spins S, i.e. S't = 3 or 5 for Co or Mn. That energy 
difference can be compared with the corresponding energy differences of isolated pairs in 
the large super-cells. Those super-cells where the magnetic ions form chains are different, 
however. Then, each magnetic ion has two nearest neighbor magnetic ions which doubles 
approximatively the previous energy difference ([2]). The exact energy difference between FM 
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and AFM states of a Heisenberg chain is slightly different, but that is unimportant for our 
present argumentation. 

The calculated exchange constants Jdd show a strong variation with U . That is illustrated 
in Fig. [T] for the Mn-doped compounds. In that case the tendency is monotonous, i.e. the 
increase of U leads to a decrease of Jdd- A similar tendency is visible for ZnA:Co with the 
exception of ZnO:Co where the LSDA exchange constant is only -1.39 meV and not of the 
order of -6. ..-7 meV like for the other compounds. This exception is due to ferromagnetic 
contributions in ZnO:Co as analyzed in Ref. [17. 



The experimental values of Jdd are known with great accuracy by magnetization step mea- 
surements or inelastic neutron scattering (see Tables HI HTj) . The comparison of experimental 
and theoretical values shows that the LSDA method strongly overestimates the exchange 
couplings. In our method the Hubbard correlation has to be chosen between 4 and 6 eV to 
obtain the correct exchange couplings. The precise value of U has also a chemical tendency. 



That was revealed in Ref. 



32 and can be explained since the compounds ZnA become less 



and less ionic in going from A=Zn to A=Te. The decrease of ionicity can be measured 
by a decrease of the charge transfer towards the magnetic ion in the series (Fig. [2]). The 
charge transfer is correlated with the calculated U value in the constrained density functional 
calculation.— Taking into account this chemical tendency we chose the U values of Tables 
m |TT]to calculate Jdd (and Na, N(3 in the next Chapter). Those values for U are shghtly 
smaller than that one calculated in Ref. 



321 since the FPLO and LMTO (linearized muffin tin 
orbitals) implementations of the LSDA+f/ method are not equivalent. Taking into account 
the restricted accuracy of our procedure we varied U in steps of 2eV. Then we obtain the 
theoretical results of Tables HI HTl which are in good agreement with the experimental values. 



IV. P-D EXCHANGE COUPLINGS 

The localized magnetic moments Sj which are provided by the magnetic ions Co^+ or 
Mn^+ interact with the spin of doped holes s. This interaction can be parametrized in the 
continuum approximation in the form: 

^ = -/5^S,s5(r-R,) , (3) 
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where the magnetic impurities are placed at R^. A similar interaction exists with the spin 
of doped electrons which is usually denoted by the parameter a. If we transform the Hamil- 
tonian into a lattice model, the interaction ([3]) becomes 

i 

with the sum over all lattice sites i which are occupied by magnetic impurities, and where 
Sj is the local spin operator of the doped hole in the lattice representation. Both parameters 
are connected by J^^ = where is the number of cations per volume (A^ = in 
the zinc blende structure). One possibility to measure A^/5 is photoemission where the hole 
in the valence band is created during the photoemission process. Another possibility is 
magnetooptics which measures the giant Zeeman effect of excitons, i.e. electron-hole pairs. 

We calculated the 'p-d exchange coupling with super-cells having impurity concentrations 
of a; = 1/4, 1/8, and 1/32 magnetic ions. The 'p-d exchange coupling leads to a valence band 
and conduction band offset between spin up and spin down IS.E'" and IS.E^ . In the case of 
weak "p-d coupling, this band offset is proportional to the impurity concentration x, i.e. it 
can be calculated in mean-field theory. That can be clearly observed in our numerical data 
and the corresponding exchange couplings are then simply given by 

^;. = iV« = -^andj;, = Ar/3 = -^, (5) 

where (S*) = {M) / {2^b) is the mean value of the local spin calculated within the ab-initio 
approach. For the Mn compounds, the calculated magnetization values are very close to 
saturation {{M) / fis = 5.00, 4.85, 4.83, and 4.85 in the series with the anions A=0, S, Se, 
and Te, respectively) but there are stronger deviations from the local value S* = 3/2 for the 
Co ones {{M)/fiB = 3.00, 2.65, 2.75, and 2.61 in the same series). 

The mean-field approach works very well for A^a which has small values in all cases. 
The reason is the small coupling between the conduction band, which is built by Zn 4s- 
4p hybridized states, with the impurity states. The calculated values are also in excellent 
agreement with the available experimental data (see Tables UTTl and HV|) . 

The situation is different for N(3. The valence band is built by the anion p-orbitals which 
have generally a large overlap with impurity states. Therefore, N(3 is much more important 
than A^a. And this tendency is increased when the lattice constant diminishes in going from 
Te to O. As a consequence, the mean-field description, and the proportionality between band- 
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FIG. 3: Valence band-edge spin splitting AE^ of ZnO:Mn (triangles) and ZnTe:Mn (squares) 
calculated with the LSDA+?7 method. The solid (dashed) lines represent the fit to the Wigner- 
Seitz model. 

offset and impurity concentration breaks down. Historically, the deviation from the mean- 
field picture was first observed experimentally for CdS:Mn.^^ In our calculations, deviations 
from the mean-field behavior are especially visible for doped ZnO and ZnS. They are mostly 
pronounced for ZnO:Mn (see Fig. [3l) where a localized state appears which means that AE"" 
tends to a constant value for a; — > (see next Section). Formally, the mean-field calculation 
of Nf3 ([5]) leads then to a divergent value which illustrates the discussed break-down in the 
most prominent way. This can also be interpreted as a crossover from the weak coupling 
to the strong coupling regime in the series going from Te to O. Since the localization is 
expected to disappear for higher impurity concentrations (visible in a band merging of the 
localized state with the valence band in the density of states) there is some justification to 
use the mean-field formula ([5]) for x = 1/4. The values of Nf3^^^ calculated in that way are 
displayed in Tables UTTl and HVl 

To resolve the deviations from the mean-field behavior a Wigner-Seitz approach was 
developed.^- We will use it to calculate NjS more accurately (see also Ref. 4^ for GaAs:Mn). 
In that theory, the valence band is described in the effective mass approximation with a spin 
dependent impurity potential. The Hamiltonian for one impurity has the form: 



^ = -2^£ + - ^;.S.s)e(6 - r) . (6) 
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Replacing the spin operator s by Sz we obtain a spin polarized scattering potential 

where a = +1(-1) =T (!) and = W - aSJ^j2 with the local spin S = 5/2 and 3/2 for 
Mn and Co, respectively. The muffin-tin radius of the scattering potential b was fixed such 
that the corresponding spheres around the cations fill in completely the space of the solid, 
i.e. (47r/3)6'^ = a^/4. The finite concentration of impurities is taken into account by the 
condition that the derivative of the wave function = vanishes at the mean radius R 

around each impurity which is determined by (47r/3)i?^ = l/{Nx). The scattering problem 
for each spin direction is easy to solve^^ and the lowest eigenvalue for spin up (down) E-^ 
{El) is given by a transcendental equation. The valence band splitting can be expressed as 

AE^ = Ei-Ei = -^Nx [Ui6{x, r/J - U^Six, v^)] , (8) 

where S{x, rja) = \Ecj/ E^^^\ is the ratio of this eigenvalue to the mean-field result 

E^^ = ^vrfe^f/^iVx . (9) 

The deviation is controlled by the dimensionless fitting parameter 



,, = ^ = 2™-[/,^-j , (10) 

where Uc is the critical potential value for the bound state creation. 

For the fit we used the gap and the band-offset calculated within the LSDA+f/ approach 
for the three concentrations mentioned above. We used the experimental values for the 
effective masses m* /m = 0.22, 0.21, 0.32, and 1.0 for ZnA, A=Te, Se, S, and O, respectively. 
These values were obtained by averaging over the transversal (light) and longitudinal (heavy) 
effective masses according to 3/m* = 2/mf + l/m;.J^ In Fig. Owe compare the weak coupling 
case (represented by ZnTe:Mn) having a linear dependence of the band-offset on the impurity 
concentration x with the strong coupling compound ZnO:Mn showing clear deviations from 
linearity. The Wigner-Seitz approach fits well our numerical data and leads to a localized 
state for ZnO:Mn. A summary of all the results is presented in Tables HIT] and HVl 



V. LOCALIZED STATE 



The Wigner-Seitz fit for ZnO:Mn results in the dimensionless coupling parameter rji = 
— 1.12 corresponding to a localized hole state. That is also directly visible in the density of 
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TABLE III: Comparison of experimental and theoretical values of the p-d exchange couplings with 
the conduction (Na) and valence {NP) bands for ZnA:Mn (A=0, S, Se, and Te). The theoretical 
results were obtained by the LSDA+t/ method, analyzed within the mean-field approximation 
{Na^^^ and Nf3^^^) and the Wigner-Seitz approach^^ (NP^^). Also given are the dimensionless 
coupling parameters r/| and r/| of the Wigner-Seitz approach. 



A^a^^P(eV) 











0.38 


-1.81 


-1.42 


0.08 


0.11 


-1.39 


-1.12 


0.07 


0.29 


-1.46 


-1.23 


0.05 


0.23 


-1.22 


-1.02 


0.04 



ZnO:Mn 
ZnS:Mn 
ZnSe:Mn 
ZnTe:Mn 



lO.l 



0.26'^ 
0.18*= 



, -2.7'' ; -3.0^ 
-1.3** 
-1.31'^ ; -1.0* 
-1.05^ ; -0.9* 



-1.12 



72 



"Magneto optical measurements, Itef 
''Photoemission spectroscopy, Ref. 
"^Photoemission spectroscopy, Ref. 
'^Magneto optical measurements, Ref. 
'^Magneto optical measurements, Ref. 
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TABLE IV: The same 


as Table [m] but for ZnA:C 


o. 








iVa^^'P(eV) 


iV/?e^P(eV) 


iVa^^(eV) 








Vl 


ZnO:Co 




1.0 (or -0.6)" ; -3.4* 


0.34 


-1.82 


-1.36 


0.34 


-0.36 


ZnS:Co 






0.21 


-2.64 


-2.24 


0.03 


-0.49 


ZnSe:Co 




-2.2^ 


0.33 


-2.50 


-1.98 


0.02 


-0.31 


ZnTe:Co 


O.Sl'' 


-3.03'' 


0.28 


-2.44 


-1.88 


0.04 


-0.34 



'^Magneto optical measure: 
^X-ray absorption, Ref. 



tnents, Ref. 47 



24 



'^Magneto optical measurements, Ref. 48 
^Magneto optical measurements, Ref. 
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states (DOS) of MnZn3i032 (see Fig. H]). A split band appears for x = 1/32, but not for 
X = 1/4. The split band indicates localization of the hole state, whereas its merging with 
the valence band for a; = 1/4 corresponds to a localization-delocalization transition with 
increasing doping. (The accurate description of this transition requires however a better 
treatment of disorder and correlation effects.) The Mn 3d majority spin states (upper part 
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of the Figure) are strongly hybridized with the valence band. Its center of gravity is located 
at about 3.5 eV below the top of the valence band. The minority Mn 3d states (lower part) 
on the contrary, are barely visible on the Figure; they start to appear at 6 eV. The split band 
is of mainly O character with a high Mn contribution. A more close analysis indicates that 
it is mainly localized on the 2p orbitals of the nearest O neighbors of the Mn impurity. As 
it is visible in the Figure, due to the isovalent impurity, the Fermi level is located just above 
the split band. Holes may be created by doping (either chemically or in the photoemission 
process). A partially filled split band corresponds to an uncompensated oxygen down spin 
which turns around the localized Mn up spin. That picture has a great analogy to the 
Zhang-Rice singlet (ZRS) state^° in cuprates. In cuprates the ZRS can qualitatively be 
described by the LSDA+t/ method in a similar manner than here. 

A localized hole state leads to several consequences. First of all, it prevents ferromag- 
netism if the doped holes are all trapped in localized states. Second, the exciton seen in 
magneto-optics is built with holes at the valence band edge and cannot be built with local- 
ized holes. However, as it is visible in Fig. HI the valence band edge is split in the opposite 
direction (apparent ferromagnetic coupling) and to a much smaller amount (about 1/3 of 
AE"). Therefore, strictly speaking, magneto-optics does not measure NjS but an apparent 
Nf3"'PP of the opposite sign and of smaller amplitude. Our LSDA+f/ calculation for ZnO:Mn 
explains this discrepancy between N(3"'Pp measured in magneto-optics and the pure antifer- 
romagnetic N(3 parameter (see Table HTB . Experimentally, the ferromagnetic sign of N(3°''^^ 
was recently unambiguously demonstrated for GaN:Fe which is not a II- VI SC, however.— 
The difference between N(3 and Nf3"'^^ can also be calculated in the Wigner-Seitz or in other 
approaches.— 

In contrast to ZnO:Mn we find no localization in ZnO:Co, but a situation quite close to 
it. In the corresponding DOS (not shown) the split band has merged with the valence band. 
It was already noted that in the LSDA calculations all 3d states are much higher in energy 
than in the LSDA+?7 (which contradicts however the photoemission measurements and is an 
artefact of LSDA). Therefore, we find hole localization in LSDA for all compounds besides 
ZnTe. Correspondingly, the \Ni3\ values are much higher (iV/3*^^ = —3.90, -2.80, -2.43, 
and -2.00 eV for ZnA:Mn with A=0, S, Se, and Te; and Np^^ = -3.86, -4.72, -4.30, and 
-4.25 eV for ZnA:Co). That contradicts the experimental data already in the weak coupling 
compounds ZnSe:Mn and ZnTe:Mn. The relevance of the LSDA+[/ approach to calculate 
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MnZn31032U = 6eV 




Energy [eV] 



FIG. 4: (Color online) LSDA+iJ density of states for MnZn3i032 (U = 6 eV, black: total DOS, 
red: partial Mn-DOS). 



NjS for ZnSeiMn was first noted in Ref. 
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which is in excellent agreement with our results. 



On the other hand, the Na values, are not very much changed by the U parameter. 



VI. DISCUSSION 



Before comparing our results with other work let us mention the limitations of our pro- 
cedure. After all, the LSDA+t/ treats correlation effects only in an approximative manner 
and neglects fluctuations. This might explain the discrepancy for ZnOiCo where state local- 
ization is very probable in view of the large difference between photoemission and magneto- 
optics.—"^ A more sophisticated method to treat correlation effects will probably refine the 
picture presented here. It means that the LSDA+f/ approach underestimates the local- 
ization tendency (and probably also the \Nf3\ values) in the strong coupling case. Other 
error sources are the limited knowledge on U, the use of the effective mass approximation 
in the Wigner-Seitz approach, which is furthermore restricted to only one valence band in 
difference to the real band structure. 

Our results show good agreement between theory and experiment for Jdd, Na and for 
Nf3 in the weak coupling regime (principally ZnTe and ZnSe). However, in the strong 
coupling case, we would like to argue that our calculated Nf] values correspond neither to the 
published ones from magneto-optics (see discussion above) nor to those from photoemission. 
Since the photoemission values of -2.7 (-3.4 eV) for ZnOiMn (ZnOiCo) were obtained in an 
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TABLE V: Photoemission data for Hubbard correlation U^jf, charge transfer energy ^eff, and 



hybridization parameter tpd, as well as the determined NP values according to Ref. 



M The p-d 



and the nearest-neig^ 
of Larson et al, Ref. 



ibor d-d exchange Jdd are determined according to the perturbation formulas 



52 





Ae//(eV) 


UeffieY) 




N/3{eY) 


Jdd(meV) 


ZnO:Mn 


7.71 


9.61 


0.80 


-2.7 


-25.29 


ZnS:Mn 


4.21 


8.41 


0.65 


-1.3 


-1.29 


ZnSe:Mn 


3.21 


8.41 


0.56 


-1.0 


-0.39 


ZnTe:Mn 


2.71 


8.41 


0.51 


-0.9 


-0.21 



indirect way using the perturbation formula of Larson et al: 



'.52 



= -ft;. 



1 



1 



(11) 



The experimental core-level photoemission spectra^^ were fitted by the configuration inter- 
action (CI) method to a MnA4 cluster (with the anions A=0, S, Se, and Te) which fixes 
the hybridization parameter tpd, the Hubbard correlation in the d shell Uejf and the effec- 
tive charge transfer energy between p- and (i-orbitals Ag// (for more details see Ref. |22I ). 
The obtained parameters are repeated in Table |V] and allow to determine Nj3 according to 
Eqn. flTT]) . The Nj3 value of -3.4 eV for ZnO:Co was obtained by an identical procedure.— 
In the same perturbation approach we may, however, also calculate the nearest-neighbor 
exchange:— 



J. 



dd 



2^2 



^pd 



1 



1 



U^ffiU^ff-A^ffY (Ae//-t/e//)3j • ^^^^ 

The calculated values are also given in Table |V] and show large discrepancies to the ex- 
perimental results (see Table [T] above) especially in the strong coupling case of ZnO:Mn. 
Similar discrepancies can be observed by determining the hybridization parameter tpd = 
{pda)/3 — 2-\/3(j5(i7r)/9 by band structure calculations.^^ These difficulties prove that the 
perturbation formulas (11111121) have a restricted applicability and have to be treated with 



mg. 



care especially for strong coup 

Being close in spirit to Ref. 120|, our results deviate nevertheless quite considerably in the 
numerical values for Nf3 which were assumed there. We found a much smaller coupling and 
we believe that the discrepancy with the published photoemission (PE) values (which are 
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about two times larger than our results) results from the non-justified use of the perturbative 
Larson formula in analyzing the PE data. As a consequence, our magnitude of the dimen- 
sionless coupling parameter 77^ = —1.12 for ZnO:Mn is much smaller than that estimated in 



Ref. 



20 



(between -2.0 and -3.3). It is highly probable, that the reduced value of \Nf3\ will 
also reduce the proposed ferromagnetic Curie temperature in ZnO:Co and ZnO:Mn provided 
that the doping level is sufficiently high to delocalize the hole states. 

The large discrepancies between different experimental and theoretical approaches for 
Nf3 in the strong coupling regime point also to the limitations of the oversimplified model 
Hamiltonian ([3]) in that limit. The p-d hybridization tpd can then no longer be regarded as 
a perturbation and the approximation of an infinite valence band width will probably lead 
to wrong conclusions. It is highly questionable that the strong coupling case can still be 
analyzed in such a manner. 

We thank Anatole Stepanov, Sergei Ryabchenko, and Roman Kuzian for useful discus- 
sions. Financial support from the "Dnipro" (14182XB) program is grateful acknowledged. 
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